Soft phonon columns on the edge of the Brillouin zone in the relaxor PbMgi/ 3 Nb 2 /303 
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We report lattice dynamical measurements, made using neutron inelastic scattering methods, 
of the relaxor perovskite PbMg 1/ / 3 Nb 2 /303 (PMN) at momentum transfers near the edge of the 
Brillouin zone. Unusual "columns" of phonon scattering that are localized in momentum, but ex- 
tended in energy, are seen at both high-symmetry points along the zone edge: Qr = { 5 , 5 , § } and 
Qm = {Ij^jO}- These columns soften at ~400K which is similar to the onset temperature of 
the zone-center diffuse scattering, indicating a competition between ferroelectric and antiferroelec- 
tric distortions. We propose a model for the atomic displacements associated with these phonon 
modes that is based on a combination of structure factors and group theoretical analysis. This 
analysis suggests that the scattering is not from tilt modes (rotational modes of oxygen octahedra) , 
but from zone-boundary optic modes that are associated with the displacement of Pb 2+ and O 2- 
ions. Whereas similar columns of scattering have been reported in metallic and (less commonly) 
molecular systems, they are unusual in insulating materials, particularly in ferroelectrics; therefore, 
the physical origin of this inelastic feature in PMN is unknown. We speculate that the underlying 
disorder contributes to this unique anomaly. 

PACS numbers: 77.80.-e, 61.10.Nz, 77.84.Dy 



I. INTRODUCTION 



The lead-based, relaxor ferroelectrics are chemically 
disordered compounds that have attracted consider- 
able scientific interest due to their remarkable struc- 
tural, dynamical, and electromechanical properties 
PbMg 1/3 Nb 2 /303 (PMN) and PbZn 1/3 Nb 2/3 3 (PZN) 
are two of the most widely studied compounds, and 
both display the broad, frequency-dependent peak in 
the dielectric response as a function of temperature that 
is characteristic of relaxors. 4 Surprisingly, even though 
both possess unusually large dielectric susceptibilities 
(e ~ 25,000), neither PMN nor PZN exhibits a well- 
defined structural transition, and thus no ferroelectric- 
ity, in zero electric field When diluted with sufficient 
amounts of PbTi0 3 (PT), which is a conventional ferro- 
electric, a long-range ferroelectric ground state is estab- 
lished and the piezoelectric coefficients grow to record- 
setting levels, making these materials prime candidates 
for industrial applications. But beyond a certain concen- 
tration of PbTi0 3 , which defines the Morphotropic Phase 
Boundary (MPB) that separates the low-concentration 
monoclinic phase from the high-concentration tetrag- 
onal phase, the piezoelectric coefficients drop precipi- 
tously ! 8 ^ 10 Whereas different ideas have been proposed 
to explain these effects, no consensus exists on the origin 



of the exceptional piezoelectric properties. 

The static properties of relaxors are highly unusual. A 
structural distortion was observed in single-crystal PZN 
using low-energy x-ray diffraction techniques^ but this 
was subsequently shown to be confined to a near-surface 
region of the crystal by another diffraction study that em- 
ployed higher-energy x-rays This localized structural 
distortion, known as the "skin effect," is enhanced by 
strong electric fields and has been confirmed using neu- 
tron strain-scanning techniques in other PZN- and PMN- 
based compounds. 13 i 14 ' 15 i 16 Despite the absence of a bulk 
ferroelectric ground state in either PMN or PZN, evi- 
dence of static, short-range, polar correlations has been 
inferred from measurements of the optic index of refrac- 
tioni 7 - and by the presence of temperature-dependent dif- 
fuse scattering that is centered on the Bragg peaksJ^ 
The refraction measurements have been interpreted in 
terms of local, nanometer-scale regions of randomly ori- 
ented ferroelectric order, often termed polar nanoregions 
(PNR), which form below the Burns temperature, Tg, 
which is nominally 620 K in PMN. 17 Various models 
of the diffuse scattering and low-energy phonons have 
been proposed; however, all models associate the dif- 
fuse scattering with the formation of polar nanoregions in 
PMN and PZN, because both appear at the temperature 
Tgff 1 ! 22 ! 23 For PT-concentrations beyond the MPB, no 
diffuse scattering is observed. As this coincides with the 
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precipitous drop in the size of the piezoelectric coeffi- 
cients, the presence of substantial elastic diffuse scatter- 
ing near the zone center is associated with the relaxor 
effects 

Whereas this diffuse scattering represents a definitive, 
static signature of relaxors, extensive research has yet 
to identify a well-defined dynamical property that is 
uniquely associated with relaxor behavior. Being chemi- 
cally disordered solids, the lattice dynamics of PMN and 
PZN are quite complex. Nevertheless, important simi- 
larities exist between the two compounds. Like the par- 
ent material PbTi03, which is a soft-mode (displacive) 
ferroelectric, both PMN and PZN exhibit phonon spec- 
tra near the zone center that are dominated by a soft 
transverse optic (TO) mode that reaches a minimum fre- 
quency on cooling and then hardens at lower tempera- 
tures J£i2yii26i27i28i29 c oncurren ^ w ith the decrease in fre- 
quency, the linewidth of this mode broadens in energy 
until it becomes heavily damped, but only for wave vec- 
tors near the zone center. This anomalous temperature- 
and wave vector-dependent damping produces a false, 
dispersion-like feature known as the "waterfall," which 
has been observed in PMN, PZN, and other relaxor com- 
pounds i 30 ' 31 Early studies speculated that the phonon 
damping resulted from a strong interaction between the 
soft mode and the PNR^2. Since then, however, the ori- 
gin of the waterfall effect has been a topic of considerable 
debate, 3 ^ 3 ^ There is now strong evidence that refutes 
the idea that the waterfall effect is caused by TO-PNR 
coupling. Low-frequency phonon studies of PMN mixed 
with 60% PbTi0 3 (PMN-60PT), a composition that is 
not a relaxor and that displays a long-range, first-order, 
cubic-tetragonal phase transition, have found a waterfall 
feature very similar to that seen in PMN.2S, Attempts 
to find a waterfall in the K(Tai_ a; Li :l ;)03 relaxor system 
have not been successful^ Thus, the heavily damped 
soft mode is not a defining characteristic of relaxors, 
and its existence does not require the presence of PNR. 
Instead, the waterfall feature probably results from an- 
other form of strong disorder common to both relaxor 
and non-relaxor ferroelectrics. A dynamic feature that 
is very probably associated with the PNR (i.e., with the 
diffuse scattering), and thus the relaxor phenomenon, is 
the strong damping of the low-energy transverse acoustic 
phonons. This damping was recently reported to have a 
strong dependence on the electric field and is absent in 
PMN-zPT compositions that lie beyond the MPB^ 

No signature of the relaxor phase has yet been iden- 
tified from studies of the zone-center soft optic mode, 
although it was recently suggested that the MPB begins 
at that PT concentration where the structural transi- 
tion temperature matches the temperature at which the 
zone-center soft phonon mode reaches a minimum. 40 The 
existence of the two temperature scales in the relaxor 
phase has been explained by random field theories i 27 ' 41 ! 42 
Finding such a feature in PMN is extremely important 
because, as is the case with soft modes in conventional 
ferroelectrics, a unique dynamical signature can provide 



information on the underlying structural distortions that 
are specific to the relaxor state. In view of the similarities 
of the zone-center soft mode dynamics throughout the 
relaxor-ferroelectric phase diagram, such as the waterfall 
effect, we believe it is necessary to search other regions of 
the Brillouin zone for a dynamic signature that is unique 
to relaxors. 



In comparison to the zone center, little experimental 
attention has been paid to the zone boundary of relaxor 
ferroelectrics, but intriguing temperature-dependent su- 
perlattice peaks have been observed in PMN using TEM 
and x-ray diffraction techniques at the reciprocal lat- 
tice vectors (|,§,0) and (~, |, — Being zone 
boundary locations these would seem unlikely to be asso- 
ciated with the relaxor transition, as any softening would 
be indicative of antiferroelectric, rather than ferroelec- 
tric, distortions. However, it has been suggested theoret- 
ically that both ferroelectric and antiferroelectric distor- 
tions can coexist, and such a concept has been applied to 
the case of SrTiC>3 in Ref. H3- The average, long-range 
structure of the cubic phase of PMN is illustrated in 
Fig. [TJ as are the high-symmetry positions of the associ- 
ated Brillouin zone. There are two high-symmetry points 
along the edge of the Brillouin zone: M = {5,5,0}, and 
R = { 5 1 5 > 5 } ■ In this paper we report the presence 
of soft phonons at both of these high-symmetry points 
in PMN. Unlike the highly damped TO phonons mea- 
sured near the zone center, the zone-boundary points dis- 
play vertical columns of phonon scattering that are es- 
pecially strong at M, and these appear at a temperature 
that coincides with the appearance of diffuse scattering 
at the zone center. Such columns of inelastic scattering 
are very unusual in insulators, and their ultimate origin 
remains unknown at this time. We focus only on the soft 
phonons in this study and do not deal directly with the 
full complexity of the elastic diffuse scattering near the 
zone boundaries reported elsewhere4£ 



This paper is divided into three sections. The first sec- 
tion provides the experiment details of the measurements 
conducted at the respective neutron scattering facilities 
located at the National Institute of Standards and Tech- 
nology Center for Neutron Research (NCNR) and Cana- 
dian Neutron Beam Centre at the Chalk River Labora- 
tories. The second section deals with the static proper- 
ties of PMN, in particular the temperature dependence 
of the elastic scattering measured at the zone boundary. 
The third section describes the inelastic measurements 
showing the presence of a temperature-dependent column 
of scattering that is associated with the temperature- 
dependent elastic diffuse scattering. We compare these 
results to other systems and, in particular, to PMN- 
60PT, which is not a relaxor, but undergoes a long-range, 
cooperative cubic-tetragonal transition, similar to that 
occuring in the parent material PbTiOs. 50 
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Pb(Mg 1/3 Nb 2/3 )0 3/ Pm3m 




FIG. 1: (Color online) a) The undistorted, average cubic unit 
cell of PMN is depicted with the Mg 1 / 3 Nb 2/ /303 octahedron in 
the center and Pb cations at the corners. The blue, red, and 
yellow spheres represent Pb 2+ , O 2- , and Mg 2+ /Nb 5+ ions, 
respectively, b) The Brillouin zone of a primitive cubic lattice 
is shown using the labels of Miller and hove— to denote the 
high-symmetry points and lines. M, R, and X represent Q — 
{I'f'O}' {I' I' I}' an d {I'O'O} reciprocal lattice vectors, 
respectively. 



II. EXPERIMENTAL DETAILS 

The measurements described in this paper are based 
on two single crystals with volumes of 3 cc and 9 cc. 
Both samples were grown using the Bridgeman tech- 
nique described elsewhere^ The large 9 cc sample is the 
same crystal that was used in a previous phonon study. 38 
Both samples have a room temperature lattice constant 
a = 4.04 A; thus one reciprocal lattice unit (lrlu) equals 
27r/a = 1.56 A" 1 . 

Neutron inelastic scattering experiments were carried 
out on the time-of-flight Disk Chopper Spectrometer 
(DCS, NCNR), the SPINS cold neutron triple-axis spec- 
trometer (NCNR) and the C5 and N5 thermal neutron 
triple-axis spectrometers (Chalk River). The experi- 
ments conducted at both the N5 and C5 thermal triple- 
axis spectrometers used a cryofurnace to vary and control 
the sample temperature between 10 K and 600 K. Exper- 
iments conducted on the DCS and SPINS instruments 
employed a high-temperature displex over a similar tem- 
perature range. To prevent any possible decomposition, 
the sample was never heated beyond 600 K and neither 
sample was ever exposed to an electric field at any tem- 



perature. Next we describe the experimental setups used 
for the various neutron scattering measurements on each 
instrument. All uncertainties (error bars) shown in the 
figures in this paper correspond to the square root of the 
measured neutron intensity. 



A. Triple-axis Spectrometers (C5, N5, and SPINS) 

Measurements on the C5 thermal triple-axis spectrom- 
eter were made with a variable-focusing (002) pyrolytic 
graphite (PG) monochromator and a flat PG (002) ana- 
lyzer. The horizontal angular beam collimation sequence 
was set to either 12'-33'-S'-29'-72' or 33'(open)-33'-S-29'- 
72' (S = Sample), and a fixed, final neutron energy 
Ef = 14.6 meV was used to define the energy transfer as 
Tiuj = Ej — Ef . A PG filter was placed after the sample to 
remove neutrons scattering from higher order reflections 
of the monochromator and a sapphire filter was placed 
before the monochromator to remove high-energy, "fast" 
neutrons. The instrumental elastic energy resolution, de- 
fined as the full width at half maximum (FWHM) of the 
elastic (fiuo — 0) peak, was ~0.9meV for both collima- 
tion sequences. The sample was aligned in the (HK0), 
(HHL), and the (H3HL) scattering planes. The (HK0) 
plane allows a complete investigation of the zone center 
acoustic phonons as described in Ref. [H. Through the 
use of the lower symmetry (H3HL) scattering plane, we 
are sensitive to the zone-boundary acoustic modes, which 
are tilt (or rotary modes) as outlined in Ref. \E^. 

Experiments conducted on N5 utilized both a flat PG 
(002) monochromator and analyzer. The final (scattered) 
neutron energy was fixed to Ef = 14.8 meV, the collima- 
tion was set at 30'(open)-40'-S'-40'-open, and a PG filter 
was placed after the sample. The instrumental elastic 
energy resolution was 1.3 meV FWHM. For these exper- 
iments the sample was aligned in the (HHL) scattering 
plane. We note that all of the inelastic data presented in 
this paper taken on triple-axis spectrometers have been 
corrected for contamination of the incident beam monitor 
as described elsewhere^ 

Measurements on SPINS used a vertically-focused PG 
(002) monochromator and a flat PG (002) analyzer. The 
final neutron energy was fixed to Ef = 4.5 meV and the 
horizontal beam collimations were set to (guide)-80'-80'- 
(open) resulting in an instrumental elastic energy res- 
olution of 0.25 meV FWHM. A liquid-nitrogen-cooled 
beryllium filter was placed after the sample to remove 
higher-order neutron contamination. 



B. Disk Chopper Spectrometer (DCS) 

Measurements on the Disk Chopper Spectrometer 
(DCS) were made at ~40 different successive crystal ori- 
entations, spaced 2° apart about a vertical axis, at 300 K 
and 600 K. The DCS instrument consists of 325 detec- 
tors in the horizontal scattering plane (or 913 includ- 
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FIG. 2: The detector trajectories, at zero energy transfer, 
for measurements on DCS with E, = 10.4 meV in the (HHL) 
scattering planes. The trajectories are plotted from ip — 10° 
to 90° with 10° spacing for clarity. The vertical dashed lines 
at H = 0.5 and 1.5 indicate the T lines along the Brillouin 
zone edge (Fig. [TJ . 



ing the out-of-plane detectors) with active dimensions 
in and normal to the scattering plane of ~3.1cm and 
40 cm, respectively. For all data presented here, only the 
central detector bank was used. The detectors are lo- 
cated 4 m from the sample at scattering angles ranging 
from 26 = 5° to 140°. Further details of the instrument 
can be found elsewhere^ The incident energy was fixed 
at 6.0 meV and 10.4 meV for measurements made in the 
(HK0) and (HHL) scattering planes, respectively, result- 
ing in instrumental elastic energy resolutions of 0.25 meV 
and 0.6 meV FWHM. The data were binned and visual- 
ized using the DCS MSLICE program^ To obtain the 
energy positions and linewidths throughout the Brillouin 
zone, we fit the data to a linear combination of simple 
damped harmonic oscillators. This lineshape is described 
in detail elsewhere (Ref. HHHil). 

The DCS experiments were conducted on the neutron 
energy-gain side, where Hlo = Ej — Ef is negative, whereas 
the bulk of our triple-axis work was conducted on the 
neutron energy-loss side where hu> is positive. The inten- 
sity of the scattering on the energy-gain side is related to 
that on the energy-loss side through detailed balanced 
Therefore, throughout this paper we retain the conven- 
tion that energy transfer is defined by Hu> = E{ — Ef and 
indicate the energy transfer to be negative for data from 
DCS and positive for the bulk of our triple-axis work. 
This experimental setup provided good energy resolution 
near the elastic position (near hu> — 0), which allowed us 
to measure the details of the low-energy acoustic and 



optic modes, and coarser energy resolution at larger neg- 
ative energy transfers, where contamination from elastic 
scattering is not a concern. 

The use of a chopper instrument may seem like an un- 
usual choice for study phonons which have dispersion in 
energy in all three directions of momentum. The use of 
DCS is different from that of a triple-axis as there are 
many more detectors whose position and scattering an- 
gle are fixed. By measuring the intensity at each detector 
as a function of time (and hence as a function of energy 
transfer), a trajectory in Q-E space is obtained. To ob- 
tain a complete map of reciprocal space as a function 
of energy transfer, many successive crystal orientations 
need to be obtained. A summary of the momentum space 
coverage at a series of crystal orentations is illustrated 
Fig. [2] The angle ip is defined as the angle between the 
[001] axis of the crystal and the incident beam. Each 
line in the figure corresponds to a particular crystal ori- 
entation and it can be seen that by rotating the sample 
orientation with respect to the incident beam, an entire 
map of Q-E space can be obtained. Throughout this pa- 
per, we have symmetrized the data around L = and 
plotted it in terms of positive values of L. 



III. ELASTIC SCATTERING 

We first examine the elastic scattering measured near 
all three high-symmetry points on the zone boundary (M, 
X, and R, as defined in Fig. Q]). Figure [3] illustrates con- 
tour plots of the elastic intensity measured with the DCS 
instrument in the (HK0) scattering plane in the vicin- 
ity of the reciprocal lattice vectors Qm = {\, f ,0) and 

Qr = (1, 1, 0). At 600K there is no observable scattering 
around M, and the T-point is characterized by a sharp 
Bragg peak. At 300 K strong diffuse scattering is seen 
near the T-point with a lineshape in q (the reduced wave 
vector, which is measured relative to the zone center in 
reciprocal lattice units) that is consistent with previous 
measurements.— A broad peak is also observed at the 
Af-point (Fig. (3Ji) indicating a second type of structural 
instability. The fact that the scattering is not resolution- 
limited in q at the M-point indicates that any static cor- 
relations are short-range and thus do not correspond to 
a well-defined structural distortion: long-range correla- 
tions would result in a resolution-limited Bragg peak like 
that shown in Fig. [3ji) . 

The temperature dependence of the zone boundary 
elastic scattering was studied on the N5 triple-axis spec- 
trometer and is illustrated in Fig.[4j The intensity at the 
X-point increases with decreasing temperature, but this 
growth matches the change in transmission through the 
crystal (measured with the straight through beam) and 
is otherwise featureless; thus it docs not represent any 
intrinsic structural change. By contrast, the intensities 
measured at both the M- and i?-points show a strong 
temperature dependence and a clear upturn near 400 K. 
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FIG. 3: Color contour plots of the elastic scattering at 300 K 
and 600 K near M (panels a) and c)) and near F (panels b) 
and d)). The data have been integrated from Twj = — lmeV 
to 1 meV and plotted in reciprocal lattice units (r.l.u.). These 
data were taken on the DCS instrument with the crystal 
aligned in the (HKO) scattering plane. 



FIG. 4: The temperature dependence of the elastic scat- 
tering measured at Q M = (§,§,0), Qr = (§,§,3), and 
Qx = (0, 0, |). Whereas a strong temperature dependence is 
observed near the M- and _R-points, little temperature depen- 
dence is observed near the X-point. These data were taken 
on the N5 triple-axis spectrometer. 



These data are consistent with the available temperature- 
dependent x-ray data of Gosula et al. (Ref. |48|), except 
that we have extended the temperature range up to 600 K 
and show that the onset is near 400 K. The enhancement 
of the elastic scattering as a function of temperature at 
the M- and i?-points is suggestive of incipient structural 
instabilities; however, it is also relatively smooth and is 
not associated with a sharp transition to a long-range, 
structurally distorted phase, which would result in a new 
Bragg peak. As the scattering is diffuse, these data infer 
the onset and subsequent growth of short-range, antifcr- 
rodistortive correlations below 400 K. This finding is sig- 
nificant because 400 K is the same temperature at which 
the zone-center TO mode reaches a minimum in energy^ 
and at which the strong, zone-center diffuse scattering, 
which is associated with the formation of static polar 
nanoregions, first appears. 24 Indeed, that the onset of 
zone center diffuse scattering occurs at 400 K in PMN 
motivated a very recent neutron scattering study of the 
Burns temperature using extremely sharp energy resolu- 
tion and resulted in a reassessed value of T^ = 420 K ± 
20 K, roughly 200 K below the value identified by Burns 
and Dacol in 1983^1 

Zone-boundary elastic scattering has been observed 
with neutrons and transmission electron spectroscopy in 
PbZro.52Tio.48O3 by Noheda et al., however, these su- 
perlattice peaks were determined to be from a minority 
domain with antiphase rotations of the oxygen octahe- 
dra about [111]. 58 X-ray studies of both PMN and PSN 
(PbSci/2Nbi/203) have also found evidence of diffuse 



scattering at the M- and i?-pointSj 48 i 59 PSN is unique 
as the perovskite B site is occupied equally by Sc 3+ and 
Nb 5+ and, therefore, ordering on this site may be ex- 
pected to induce at least short-range correlations at ei- 
ther the M- or i?-points. The -| : | stoichiometry on 
the B site of PMN does not allow such ordering, which 
means that short-range correlations at either Mori? are 
not necessarily expected. Consequently, it is important 
to search for any phonon anomaly at the M- and i?-points 
to determine the symmetry of the distortion, as well as 
to determine whether the elastic scattering originates in 
the bulk, and thus is an intrinsic effect, or whether it is 
limited to the surface or "skin" region and is influenced 
by strain or local compositional inhomogeneities. 



IV. INELASTIC SCATTERING 

Phonon dispersion measured below 30 meV, for both 
transverse and longitudinal polarizations, are plotted in 
Fig. [5] These data are based on the experiments de- 
scribed in the previous section and represent a compila- 
tion of triple-axis and time-of-fiight measurements. Pre- 
vious experiments have focused on phonon branches that 
include the T-point as well as on the temperature depen- 
dence of the soft transverse optic and acoustic modes. 
The data for those measurements are summarized at 
room temperature in panels a), b), and e) for modes 
propagating along the [111], [110], and [001] directions, 
respectively. Whereas we will not discuss the longitudi- 
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FIG. 5: A compilation of the room temperature (T = 300 K) 
phonon dispersion curves in PMN measured along high sym- 
metry directions within the Brillouin zone. The locations of 
the incipient soft modes are denoted by grey columns: M 
(zone boundary), Y (zone center), and 7? (zone boundary). 
These data were obtained from experiments conducted on 
the DCS chopper instrument, and the N5, C5, and SPINS 
triple-axis spectrometers. Panels a)-e) display the transverse 
phonon dispersions and panels f)-h) show the longitudinal 
dispersions. The solid lines are guides to the eye. 



nal scattering in detail in this paper, we summarize our 
current data set to define the energy scale of the phonon 
scattering in PMN in panels f), g), and h): the energy 
scale is similar to that reported for PbTiC>3. 

The dispersion curves reveal several possible incipi- 
ent soft mode positions, these being where the optic 
modes reach a minimum energy as a function of reduced 
wave vector. Aside from the soft TO mode at the T- 
point, which has been studied in detail by several groups, 
optic mode minima are seen at the M- and i?-points; 
these minima correspond to the locations where strong, 
temperature-dependent elastic scattering was observed, 
as described in the previous section. In this paper we 
concentrate on the phonons and the dispersion curves 
along the T-line, which is the zone-boundary edge that 
links the high-symmetry M- and i?-points (see Fig. [5] 
c)). We next describe the temperature and wave vector 
dependence of the soft modes, and we assign these to spe- 
cific optic modes based on our measurements conducted 
with time-of-flight and triple-axis spectrometers. More 
importantly, we show that these soft modes cannot be 
the tilt instabilities in cubic perovskites, which are also 
associated with phonon instabilities at these same points 
along the edge of the Brillouin zone. 
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FIG. 6: Inelastic scattering measured at 300 K and 600 K 
along the T-edge of the Brillouin zone (Qt = |, |,L), where 
L = 1 is an A/-point, and L = ^, | are i?-points, indicated 
by the vertical dashed lines. These data are integrated over 
(<5H, <5H, 0) = (± 0.05, ± 0.05, 0) in r.l.u. and were taken on 
the DCS instrument with E,:=10.4meV. The negative energy 
transfer on the y-axis corresponds to an energy gain of the 
neutron as explained in the text. The intensity has been nor- 
malized by the Bose factor to put the intensities on an equal 
scale. 



A. DCS Time-of-Flight Spectra 

We summarize the unusual zone-boundary phonon 
scattering in Fig. [SJ which displays DCS spectra recon- 
structed along the T-line (§,§,L) with Qm = (§ 5 §>1) 

and Q R = (§,§,3) and (§,§>§) at 600K and 300K - 
These color contour plots illustrate three key features re- 
garding the nature of the lattice dynamics near the zone 
boundaries. 

First, a soft mode is evident at the M-point in PMN 
given the substantial increase in spectral weight that oc- 
curs at low energies upon cooling from 600 K to 300 K, 
which results in the column of scattering shown at Qm = 
(§, §, 1) in Fig.H At 600 K (panel a)) the data show the 
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presence of a broad mode centered near 4 meV. At 300 K 
(panel b)) this mode softens considerably on cooling, ex- 
tending all the way to the elastic position, but it does so 
in such a way that it essentially becomes a continuum of 
scattering over energy transfers from 0-5 meV. 

Second, another soft phonon is observed at the i?-point 
Qr = (I, |, |). At 300 K there is evidence of a column, 
but it is far weaker and much narrower in momentum 
than that at Qm = (|,|,1). This column appears to 
extend down to the elastic position clearly at 300 K. It 
is also important to note the absence of any column at 
Qr = (§> |) §)■ We attribute this to the phonon struc- 
ture factor, and we use this observation later when as- 
signing the symmetries to the incipient soft modes at M 
and R. 

A third feature, which is common to both the M- 
and i?-points, is that the strong softening of these zone- 
boundary modes has an unusual dispersion that is local- 
ized in momentum but broad in energy, which is what 
gives rise to the appearance of a column of scattering 
in Q-E space. This feature is anomalous amongst the 
ferroelectrics and perovskites, where soft modes are usu- 
ally described by a well-defined dispersion that is charac- 
terized by a single, temperature-dependent value of en- 
ergy for a given momentum transfer, as seen near the 
r-point in several studies of PbTi0 3 , PMN and PMN- 
60PT. Since the soft modes observed at both the Al- 
and i?-points exhibit a continuum of energies at these 
specific momentum transfers, they cannot be understood 
by such a well-defined function. 
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FIG. 7: Constant-Q scans measured at the Qm = (|,|,0) 
M-point using triple-axis spectrometers. Panels a) and c) 
show scans at 600 K and 300 K, respectively, that demonstrate 
a dramatic softening on cooling. Panel b) is a color con- 
tour map, compiled from many constant-Q scans measured 
on SPINS, that illustrates the shift of intensity between high 
and low energies as a function of temperature. The data were 
taken in the (HK0) scattering plane. 



B. Triple- Axis Measurements at the M- and 
i?-Points 

To determine the origin of the inelastic scattering, un- 
derstand the temperature dependence, and assign the 
soft mode to a particular symmetry, we have conducted 
a series of triple-axis experiments in the (HK0), (HHL), 
and (H3HL) scattering planes, each of which provides 
access to different types of phonon modes. 

We first examine the M-point phonons. Fig. [7] displays 
two constant-Q scans in which the scattered intensity was 
measured at a fixed momentum transfer Qm=(^,§,0) 
while scanning the energy transfer E = hui = Ej — Ef at 
600 K (panel a)) and 300 K (panel c)). Whereas two well- 
defined modes are present at 600 K, the 300 K data show 
a broad distribution of scattering resembling a contin- 
uum that extends up to ^7meV. The lineshape used to 
fit these data is based on the sum of two damped, simple 
harmonic oscillators and has been described elsewhere*^ 
Panel b) illustrates the evolution of the phonon scattering 
with temperature at Qm — (^, §, 0) through a color con- 
tour plot compiled from constant-Q scans measured from 
10 K to 600 K. At high temperature little scattering is 
seen below ~4meV; at lower temperature (near 300 K) 
a clear continuum of scattering develops with intensity 
extending from 4 meV to the lowest energy transfer in- 



vestigated. The large intensity below Tioj — 0.5 meV is 
attributed to pure elastic scattering that is sampled due 
to the energy resolution of 0.25 meV. For temperatures 
below 300 K, a decrease in the scattering is observed in 
the energy range 0-4 meV. This behavior is consistent 
with that of the soft mode observed near the F-point 
and illustrates the transfer of spectral weight from the 
inelastic channel to the elastic channel, discussed in the 
previous section. 

The temperature dependence of the M-point inelastic 
scattering was characterized using constant-energy scans. 
This method was chosen over the conventional constant - 
Q method because the column is a spectral feature that is 
well defined in momentum, but broad in energy. There- 
fore, a scan of momentum transfer Q, at a fixed energy 
E, will result in a well-defined peak centered at the zone 
boundary. 

Fig. [8] illustrates constant energy scans measured at 
Tiuj = 1.55 meV over 10 K to 600 K. The data were 
taken on the C5 triple-axis spectrometer using a fixed fi- 
nal neutron energy Ef = 14.6 meV. Panel a) summarizes 
the peak intensity, taken as the amplitude of a Gaus- 
sian fit to the correlated scattering in momentum across 
the M-point, as a function of temperature. On cooling, 
the intensity at 1.55 meV increases, peaks at T ~ 400 K, 
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FIG. 8: A summary of constant-E scans measured at an en- 
ergy hio = 1.55 meV across the Q.m=(~,|,0) M-point col- 
umn. Panel a) plots the intensity of the peak observed in 
these scans as a function of temperature, and representative 
scans are displayed at 10 K, 350 K, and 550 K in panels b), 
c), and d). The data were taken on the C5 (Chalk River) 
thermal triple-axis spectrometer using E/ = 14.6 meV in the 
(H3HL) scattering plane. 



FIG. 9: A summary of constant-E scans conducted at Tujj = 
3.1 meV through the .R-point column at Q_r=(|, |, |). Panel 
a) plots the peak intensity observed in scans through the in- 
elastic column, and representative scans are displayed at 10 K, 
350 K, and 600 K in panels b), c), and d). The data were taken 
on the N5 (Chalk River) thermal triple-axis spectrometer with 
fixed Ef — 14.8 meV in the (HHL) scattering plane. 



and rapidly declines, becoming nearly unobservable by 
10 K. Examples of the fits at several temperatures are 
illustrated in panels b)-d). 

Whereas our preceding discussion has focused around 
the M-point, the Q-E data taken on DCS reveal a sec- 
ond column at the i?-point. Fig.[9]shows constant-energy 
scans at Tuo — 3.1 meV. Panel a) shows the temperature 
dependence of the column. Panels b)-d) display rep- 
resentative constant-energy scans fit to Gaussian peaks 
centered at R. The scattering at R (Fig. [9} shows a tem- 
perature dependence on cooling nearly identical to that 
observed at M (Fig. [5]), reaching a maximum intensity 
at 400 K and declining at lower temperatures. The tem- 
perature dependence of the scattering observed at both 
R and M cannot be accounted for by the Bose factor 
alone; if this were the case, then the intensity at a fixed 
energy transfer would continue to grow with increasing 
temperature. Therefore, our results show that there are 
incipient soft modes at both M and R that reach a min- 
imum frequency (hence giving maximum intensity along 
the energy column) at ~400K. 

It is important to work out whether the soft modes at 
M and R are acoustic (tilts of the octahedra) or optic 
modes. It is difficult to answer this question using scans 



exactly at these high-symmetry points, where the soft 
modes form a broad continuum in energy. We therefore 
found it useful to conduct constant-Q scans at wave vec- 
tors slightly displaced from M and R, where the acoustic 
and optic modes could be resolved at all temperatures. 

Constant-Q scans near the i?-point are presented in 
Fig. [I0l where the solid line is the result of a fit to two 
damped harmonic oscillators. Two distinct modes are 
observed, with the lower-energy mode showing little tem- 
perature dependence while the higher-energy mode hard- 
ens upon cooling to 100 K. As this upper mode is temper- 
ature dependent, we associate the columns of scattering 
at the i?-point with this optic mode, and not the tilts. 
This point is confirmed directly using structure factors 
and group theory analysis in the final section of this pa- 
per. We have not pursued the temperature dependence 
of the frequencies at high temperatures (to observe the 
softening) due to the possibility of sample breakdown. 
This softening would require measurements over a very 
broad temperature range as the rate of softening with 
temperature is predicted, and has been observed in con- 
ventional ferroelectrics, to be much less than the rate of 
recovery at low temperatures. 

Similar scans were also conducted near M and are pre- 
sented in Fig. QT] The solid line in the figures is a fit 
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FIG. 10: Constant-Q scans conducted at a wave vector dis- 
placed from the i?-point at 600 K, 300 K, and 100 K. The 
data illustrate a softening of the optic mode that recovers to 
an energy of ~15meV at 100 K. The data were taken at the 
C5 (Chalk River) thermal triple-axis spectrometer with fixed 
Ef = 14.6 meV. 



to the sum of two damped harmonic oscillators. For 
scans at Q = (i, §, 0.125), which lies close to M, two 
low-energy phonon modes could be clearly observed at 
600 K (panel a) of Fig. [□}. At 300 K (panel b)) only 
a continuum of scattering is observed, but a recovery 
of a higher energy branch is apparent at 100 K (panel 
c)), where the higher energy mode recovers to an energy 
transfer of Ticu = lOmeV. Whereas the interpretation 
of the low-energy scattering in these scans is ambiguous, 
because the M-column is intense and distributed in mo- 
mentum (see FigEJ) and there are several modes that may 
be close in energy (see Ref . [6fj| ) , these scans do show that 
there exists a higher-energy mode that is temperature 
dependent. This implies that an optic mode contributes 
strongly to the column of scattering at this wave vec- 
tor. We rule out the contribution from zone-boundary 
acoustic, or tilt, modes through a discussion of structure 
factors and group theory presented later in the paper. 




FIG. 11: Constant-Q scans conducted at a wave vector dis- 
placed from the M-point at 600 K, 300 K, and 100 K. The 
data illustrate a softening of the optic mode, which recovers 
to an energy of ~10meV at 100 K. The data were taken at 
the C5 (Chalk River) thermal triple-axis spectrometer with 
fixed E f = 14.6 meV. 



C. Comparison with a "Conventional" 
Ferroelectric (PMN-60PT) 

We now make a brief comparison with PMN-60PT, 
a conventional ferroelectric that undergoes a first-order, 
cubic-tetragonal transition at 540 K. The large con- 
centration of PT places PMN-60PT beyond the MPB, 
which divides relaxor ferroelectrics, defined by broad 
temperature- and frequency-dependent peaks in the di- 
electric susceptibility, from conventional ferroelectrics, 
defined by sharp peaks in the dielectric susceptibility 
and well-defined structural transitions. The transition 
in PMN-60PT is driven by a soft TO mode and has been 
studied as a function of temperature. The soft TO mode 
is the same mode that softens in relaxors, such as PMN 
and PZN, and that possibly gives rise to the static diffuse 
scattering. 

A summary of the phonon dispersion, measured in 
the (HK0) scattering plane, is illustrated in Fig. for 
T = 540 K (panels a)-c)) and T = 800 K (panels d)- 
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FIG. 12: Constant-Q scans measured near the M-point in 
PMN-60PT at 800 K and at 540 K, in the vicinity of the cubic- 
tetragonal structural transition. The data were taken on the 
N5 (Chalk River) spectrometer with fixed E/ = 14.8 meV. 



V. GROUP ANALYSIS 

Having identified that the columns at R and M are 
the result of soft optic (not acoustic or tilt) modes, we 
now investigate which atoms contribute to these modes 
through the use of group theory and structure factors. 
The history of the lattice dynamical studies of perovskites 
is so long that various systems have been used to label 
the eigenmodes of the phonons, including a system used 
by Cowley^ Miller and Love^ and Bouckaert, Smolu- 
chowski, and Wigner (BSW)^ The translation between 
the labeling system of Cowley and BSW is given by 
Boyer.— The established setting of the crystal in most 
lattice dynamical studies of perovskites is one where the 
octahedron sits in the center of the unit cell, centered 
on the Wyckoff &-sites with O on the c-sites, and Pb on 
the corner a-sites (Fig. Q]a). This contrasts with most 
discussions of the structure and group theory underlying 
the tilt transitions in perovskites, in which, the octahe- 
dron conventionally sits on the corner of the unit celh 74 i 75 
This difference in the cell setting causes the mode la- 
bels, even within the same labeling system, to change. 
In order to make a comparison to the majority of the 
lattice dynamics studies we use the former setting, typ- 
ically used for dynamical studies. The decompositions 
were performed with SMODES^ and the compatibil- 
ity relations were derived with the Isotropj*^ code. We 
calculated the structure factors based on the mode eigen- 
vectors derived from Isotropy^ and visualized them with 
Isodisplace^ We use the labels of Miller and Love4^ 



f)) near M. At the zone boundary (panels c) and f)), 
a well-defined phonon mode is observed near 6meV at 
both 800 K and 540 K. The phonon linewidth is smaller 
than the energy of the phonon and therefore the peaks 
represent underdamped propagating phonon modes. The 
phonon energy is unchanged within error on cooling from 
800 K to 540 K, a temperature very near to the structural 
transition temperature. The phonon peaks are well de- 
fined and no column or continuum of scattering is ob- 
served near M (q = 0.5). Constant-Q scans at q = 0.3 
(panels a) and d)) and at q = 0.4 (panels b) and e)) show 
that, in contrast to PMN, well-defined phonon modes are 
observed. No continuum scattering has been reported in 
pure PbTi0 3 (PT) and we therefore assume it to be ab- 
sent as it is in PMN-60PT^i 

It seems likely that the soft columns at the high- 
symmetry points of the Brillouin zone edge are only 
present in the relaxor phase and possibly only for PT 
concentrations below the MPB. This will need to be ver- 
ified through a study of several other lead-based relaxor 
ferroelectrics, particularly near the MPB concentration, 
which for PMN-xPT is x = 32%. 



A. Mode Decomposition and Structure Factor 
Analysis 

The Brillouin zone is shown in Fig. [TJa, along with la- 
bels of high-symmetry points and lines, using the sys- 
tem of Miller and Love42 The total decomposition of 
the atomic motions into irreducible representations was 
calculated at T, M, and R and compared to that de- 
rived by Cowley^S, A compatibility diagram is shown in 
Fig. Q21 which also relates our notation to those listed 
previously. To remove ambiguity as to which modes 
at M are linked to those at R, published, measured, 
and calculated dispersion curves of SrTi0 3 ,^ 2 RbCaF 3 j2* 
CaSi0 3 ^ KNb0 3f Sl PbTi0 3 ^>££ BaTi0 3) ^ KTa0 3r ^ 
KMnF^ and PbZrOg 7 ^ were consulted. There are rel- 
atively few i?-modes. A common feature of the low- 
frequency lattice dynamics of all the perovskites listed 
above is that the two lowest frequency i?-modes are R§ 
and iij , whereas the highest frequency mode is i?^ . Ow- 
ing to the number of modes at M, there are fewer general 
observations that can be made about the spectra at this 
wave vector. 

The i?-point is the corner of the Brillouin zone, which, 
being a cube, lies at the convergence of three T zone edges 
(Fig. [TJd). Therefore, in general, the i?-point irreducible 
representations are of higher dimensionality than those at 
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FIG. 13: The compatibility ol the irreducible represen- 
tations of phonon modes propagating along M-T-R is il- 
lustrated. The two lowest modes (R^ and R^) are the 
lowest energy modes in SrTiOs,— RbCaF 3l — CaSiOa,— 
KNbOs,— PbTiOa,— & BaTiOs,— KTaOs,— , and PbZr0 3 .— 
ML stands for the labeling system of Miller and Love» 9 , C for 
that of Cowley^ and BSW that of Bouckaert, Smoluchowski, 
and Wigner.— The symmetries of the zone boundary columns 
are highlighted. 



M; on moving from R to T the loss of degeneracy means 
that, with the sole exception of the highest-frequency i?^ 
mode, each i?-mode splits into two T-branches terminat- 
ing in a single M-mode. Fig. [Tfl plots contours of phonon 
intensity of the column near the i?-point at 300 K, based 
on data taken on the DCS instrument. The rr-axis is 
the Q = (|,|,L) direction, i.e., the T-line denoted in 

Fig. Hb. The plot shows the soft mode at Qr = (|, |, |) 
splitting into the two T-branchcs: one low-energy mode 
that resides between 5-8 meV, and another mode that 
rises rapidly in energy to over 20meV. 

One famous feature of the low-frequency dynamics of 
perovskites is the tilt modest All of the commensurate 
tilt structures of perovskites can be related to combina- 
tions of various dimensions of two zone-boundary acous- 
tic modes at M and i? i 74 i 75 The T-line that joins these 
points has also been shown to have tilt modes^ and the 
whole M-T-R line has been observed to soften in some 
cases such as SrTiOai 7 -^ The tilt mode at R is RJ , which 
splits into one low-frequency Ti tilt mode, terminating in 
the tilt mode, and one non-tilt T5 mode, which rises 
rapidly in frequency, terminating in a mode of symme- 
try M 5 ~ (Fig. [T3")) . Similarly, the low- frequency ii7 optic 
mode splits into the Tg- and Ti-branches, which termi- 
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PMN, T=300 K, DCS 
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FIG. 14: Intensity contours of the phonon scattering near 
the f?-point measured at T = 300 K, plotted on a logarith- 
mic scale. A higher-energy mode originating from the 7?- 
point (labelled T5) is clearly seen extending up to at least 
-~20meV energy transfer. The data were taken on DCS with 
Ei=10.4meV. 



nate in and M% , respectively: there is, therefore, the 
potential for an anti-crossing (Fig. [T3|) of the T^-branches 
in the low-frequency dynamics, as is observed near R in 
RbCaF^ and CaSi0 3 ^ These represent the likeliest 
candidates to be observed in the lowest frequency part 
of the spectra. To avoid relying solely on these general 
observations, we have examined the structure factors of 
all the M and R modes so as not to exclude other possi- 
bilities. 

The inelastic structure factor is given by the following 
formula for the m-th mode, 



M3) = E(3-^f J&»e*K (1) 

where Q is the total momentum transfer, ^ m ,n is the 
eigenvector for the m-th eigenmode, and r n is the po- 
sition, b n the scattering length, and m n the mass of the 
n-th atom. In this equation n sums over all atoms in 
the unit cellj 7 -^ Since we are working at a fixed tempera- 
ture and for comparably small momentum transfers, we 
have ignored the Debye- Waller factor. Using this for- 
mula and the eigenvectors obtained from Isotropy^ we 
find that tilt distortions at Qr = (§, §, 5), where we ob- 
serve the strong column of scattering, have identically 
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Eigenvectors, R point columns 
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FIG. 15: The possible eigenmodes resulting in the _R-point 
zone boundary column. The blue, red, and yellow atoms cor- 
respond to Pb 2+ , 2 ~, and Mg 2+ /Nb 5+ ions, respectively. 
Panel a) illustrates the tilt displacements about the x-axis 
corresponding to a single component of all neighbor- 

ing octahedra above and below would be tilted out-of-phase. 
Panel b) illustrates the displacements of a single-fc distortion 
of the 7?J optic mode. Based on a comparison with SrTi03 
and a structure factor analysis, we suggest that the column 
is associated with displacements of the form shown in b), i.e., 
R 4 . 



zero structure factor in the zones in which the DCS spec- 
tra were measured. Therefore, the modes contributing to 
the columns must instead be zone-boundary optics, i.e., 
dynamic antiferroelectric distortions. 

By comparison with many other cubic perovskites, the 
likeliest next-lowest frequency optic mode at R is 
(Fig. [T5]) . This is particularly clear in the calculations 
of Garcia and Vanderbillj^ for PbTiC>3 that have been 
confirmed experimentally by Tomeno et al§^ In addition 
to the general observation about the order in frequency 
of modes at R, calculations of the structure factors of 
i?-point phonons show that a single-A; structure of i?7, 
in which Pb and O are displaced, can lead to near-zero 
neutron structure factors at those i?-points for which all 
indices are equal, as observed in our scattering data at 
ll 11) (3 3 3) t 

V2 ' 2 ' 2'' 2 ' 2' 2' 1 

Having assigned the i?-column based on a combination 
of group theory and structure factor analysis, and hav- 
ing eliminated the tilt modes from consideration, we now 



assign the displacements associated with the M-column. 
As illustrated in Fig. QTJ a single T-branch appears to be 
interrupted by columns at the high-symmetry M- and 
i?-points. As the Ts-branch, due to anti-crossing, typ- 
ically rises rapidly in energy in other perovskites, the 
symmetry of the M-column is probably compatible with 
the Ti-branch. Assuming this, we assign the symmetry 
of the M-column phonon to Mg~ . This conclusion is cor- 
roborated by structure factor calculations. We observe 
a soft phonon column at M-points of all forms, includ- 
ing (|, |,0) and (i, |,0). The inelastic structure factor 
of the M^~ phonon is finite at all M positions and grows 
with increasing \Q\. The eigenvectors of Mg~ consist only 
of displacements of Pb ions along the (100) directions. 
We note that this conclusion is in agreement with calcu- 
lations presented in Ref. [79|, which predict a low energy 
phonon mode consisting of antiphase displacements of 
Pb. 



VI. DISCUSSION 

It is interesting to compare the lattice dynamics of re- 
laxors to those of conventional ferroelectrics. Two no- 
table features are present in PMN. The first is the 
"waterfall" effect, which is characterized by a broaden- 
ing in energy of the low-energy TO mode for wave vec- 
tors near the zone center. This phonon anomaly is also 
present in PMN-60PT, which exhibits a sharp structural 
phase transition, and in rclaxor ferroelectrics having PT 
concentrations below the critical concentration of 32%, 
which defines the MPB. As the exceptional dielectric 
properties of relaxor ferroelectrics are absent in PMN- 
60PT, the "waterfall" effect cannot be linked to relaxor 
behavior. The second notable feature is the presence of 
columns of inelastic scattering observed at the M- and 
i?-point zone boundaries in PMN. These columns exhibit 
a similar softening near 400 K, where strong diffuse scat- 
tering is onset (as shown in other experiment s 24 ! 57 ), but 
they are not present in PMN-60PT. Therefore, whereas 
the zone-boundary soft phonons may not be the direct 
cause of the unique dielectric properties in PMN, they 
may be indirectly coupled and indicative of the underly- 
ing relaxor ferroelectricity. 

There is no obvious physical reason for these columns 
of scattering. Similar columns have been observed in 
several other systems; examples include TTF-TCNQ 
(Ref. M), K 2 Pt[CN) 4 Bro.3 ■ 2.3 D 2 (Ref. MM), 
K0.3M0O3 (Ref. [83), and body-centered cubic zirconium 
(Ref. Hi]). In the above cases, the near vertical rod 
of phonon scattering has been attributed to a Kohn 
anomaly or a coupling between electronic and lattice de- 
grees of freedom^ 5 - Columns of scattering in insulators, 
where no such coupling can exist , are much rarer but have 
been observed in CaC0 3 (Ref. M) and KCN (Ref. 13). 
Both of these systems display molecular rotations and the 
transitions are order-disorder in nature; here the columns 
were interpreted in terms of a coupling between molec- 
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ular relaxational modes and harmonic lattice vibrations. 
Theories describing such coupling have been proposed 
by Michel and Naudts (Ref. l88l ) and similar models have 
been used by Shapiro et al. (Ref. [89t ) to describe the 
central peak in SrTiC>3. Therefore, columns of scatter- 
ing have been associated with a coupling of the lattice 
(harmonic degrees of freedom) with another degree of 
freedom (either electronic or molecular for the examples 
cited) . 

The soft phonon columns in PMN are therefore very 
unusual. The lineshapc, which is localized in momen- 
tum but broad in energy, cannot be associated with an 
electronic process (such as a Kohn anomaly) as PMN 
is an insulator and therefore there is no Fermi surface. 
Nor can they be directly associated with a coupling to 
molecular rotations, as no such local molecular ions exist 
in the PMN structure. One possibility is of a local rat- 
tling of atoms, which is directly related to the disorder in 
PMN. The local disorder and random fields caused by the 
charge difference between Nb 5+ and Mg 2+ may introduce 
a number of local states that can be occupied (in analogy 
to molecular solids) , to which the lone pair displacement 
of Pb 2+ ions can couple. Such structures have been pro- 
posed to understand powder diffraction measurements in 
PZN (Ref. HI and PMN (Ref. i) as well as NMR mea- 
surements (Ref. [9Tb . The role of Pb 2+ ion displacements 
in the unusual piezoelectric properties was also proposed 
based on PDF analysis^ Therefore, as the same atoms 
can be on several sites in different unit cells, this may act 
as a means for a relaxational mode that can then cou- 
ple to the harmonic modes at high energies. Columns 
are not observed in PMN-60PT, which lacks any tem- 
perature dependent zone-center or zone-boundary elastic 
diffuse scattering, but are present in PMN, where strong 
zone-boundary elastic diffuse scattering is present, infer- 
ring that static, short-range correlations are required for 
columns to form. The columns of scattering seem to rep- 
resent a series of nearly degenerate states, which must 
be related to the formation of a non-zero polarization, 
as evidenced by the soft mode at the T-point. The pres- 
ence of several nearly degenerate states is also evident 
from the complex structural phase diagrams and proper- 
ties observed as a function of PT-doping.— We note that 
this idea is completely heuristic, and thus it will be inter- 
esting if future calculations and theories can reproduce 
these columns. 

The temperature dependence of the soft phonon 
columns and the fact that the softening occurs at nearly 
the same temperature as the minimum in the frequency 
of the zone-center TO mode and the onset of the diffuse 
scattering is strongly suggestive that the zone-boundary 
columns maybe a signature of the relaxor state. This is 
confirmed by the absence of of these columns in PMN- 
60PT, a well-defined ferroelectric. The only conclu- 



sive means of proving the direct connection with the 
relaxor state is a detailed study as a function of PT- 
concentration. 

Our results point to an underlying antiferroelectric dis- 
tortion in the relaxor PMN that competes with the fer- 
roelectric distortion, characterized by strong, zone-center 
diffuse scattering and a soft optic mode. The presence 
of strong random fields and competing interactions pro- 
vides an underlying physical reason for the glassy nature 
observed using dielectric measurements. How these com- 
peting interactions change with doping and, in particular, 
near the MPB will be a topic of interesting future study. 



VII. CONCLUSION 

We have discovered the existence of very unusual soft 
modes in PMN at the Qr = {5,^,5} and Qm = 
{i,i,0} zone boundaries. The phonon modes have an 
anomalous dispersion, arc localized in momentum, but 
broad in energy resulting in a column of scattering. The 
columns of scatter is present in PMN (a relaxor ferroelec- 
tric) are absent in PMN-60PT (a well-defined ferroelec- 
tric). The physical origin of the unusual dispersion is not 
understood and cannot be accounted for by a coupling 
between electronic and lattice degrees of freedom, nor the 
presence of a molecular relaxational mode. The preva- 
lence of many inequivalent local atomic sites may result 
in a relaxational mode, in analogy to order-disorder sys- 
tems. It is possible that a strong coupling between zone- 
center and zone-boundary distortions may be a general 
component of relaxor ferroelectrics, which requires ver- 
ification for further compositions on either side of the 
MPB. 
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